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In this Letter we consider the adiabatic charge transport through a normal mesoscopic sample 
sandwiched by superconductors without modulation of local chemical potentials. The deformation 
of coherent quasiparticles in the normal metal in the presence of periodically changing phases in the 
superconductors leads to a charge transport. Both the magnitude and the phase dependence of the 
charge transferred per period through the sample strongly depend on temperature. 



Recently adiabatic charge transport in mesoscopic 
samples generates much interest both experimentally and 
theoretically . In the presence of slowly and period- 
ically changing external fields, the spatial electron den- 
sity profile follows periodical modulation of local chem- 
ical potentials. In the presence of two a.c. charging 
gates, {gi{t) , g2{t)} , this leads to a d.c. adiabatic charge 
transport across the sample in the classical limit. The 
amount of charge transport equals the amount of "flux" 
of a topological field ttu (g) threading the area enclosed 
by the trajectory C swept by g = {gi{t), g2it)} in a 2 
dimensional parameter space A^2- At low temperature 
when the dephasing length is longer than the sample size, 
electrons remain coherent. 7ri2(g), unlike in the classical 
limit where it is a constant field, develops very rich struc- 
tures and is a sample specific, random function of g in 
M2 space. Due to quantum interference effects, the adi- 
abatic charge transport at low temperature therefore is 
of a random sign and mag nitude 

It still remains a question how the many body correla- 
tion enriches the structure of the topological field 7ri2 (g) 
and affects the adiabatic charge transport. In this let- 
ter, we exam an example where a mesoscopic sample is 
in contact with superconductors. At T ^ A(the super- 
conductor energy gap), electron-like excitations in these 
samples are reflected into hole-like excitations through 
Andreev reflections at superconductor-metal(SN) inter- 
faces. Quasiparticles inside the normal metal at low tem- 
peratures are superpositions of electrons and holes, with 
zero effective charges [^-^ . For a superconductor- normal 
metal-superconductor (SNS) junction, the resultant ex- 
citation spectrum in the metal has little resemblance to 
that of a normal Fermi liquid: An energy gap of many 
body nature is opened up at Fermi energy, which remains 
finite while the single particle level spacing goes to zero 
[^-||. The transport coefficients are also renormalized. 

When a voltage is applied across the junction, the 
phase difference is changing according to Josephson re- 
lationship. For the geometry shown in Fig.l, the phase 
differences of two SNS junctions are given as 

Xl.2W-^+<5l,2 (1) 

with di,2 as the phase differences when V — 0. Electrons 
experiencing Andreev reflections would acquire different 



phases at SN interfaces. Thus, electron-hole superposi- 
tions of quasiparticles depend on phase differences of su- 
perconductors. The energy gap in the quasiparticle spec- 
trum decreases from a finite value to zero when the phase 
difference between two superconductors is changing from 
zero to TT |^-|| . Quasiparticles of energy e at time t can 
be represented by {vi{e, t), V2{e, t)), with t'oi^i(e, t) as the 
average density of states in region 1, voU2{e,t) in region 
2. (The density of states in the normal metal outside 
region 1,2 equals its bulk value i/q, doesn't depend on 
Xi,2it)-) As time evolves, xi,2 changes and i^i_2(e,i) os- 
cillates with a frequency 2eV/h. When Xi(2) is far away 
from from tt, quasiparticles far below the energy gap ex- 
tend only outside region 1(2). When Xi(2) approaches tt, 
the gap diminishes; the low energy quasiparticles start 
to extend into region 1(2). This leads to a continuous 
deformation of quasiparticles of energy e, characterized 
by the motion of (i^i(e, t), ;/2(e, t))(see Fig. 2). It is worth 
emphasising that unlike the situations in [^]J^,^, the lo- 
cal chemical potential in this case remains constant as 
a function of time and the deformation of many body 
wave functions is possible only when the phase rigidity 
is present |Q. For this reason, we call the phenomena 
discussed here as phase pumping. 

In the following, we show that deformation of coherent 
quasiparticle wave packets via phase differences eventu- 
ally leads to a d.c. current in this case. The adiabatic 
charge transport in the presence of the changing phase 
differences of superconductors is determined by the sensi- 
tivities of the diffusion constants and the density of states 
to boundary conditions. At low temperatures, the topo- 
logical field is localized in the parameter space, leading 
to a very singular phase dependence of d.c. charge trans- 
port. The charge transport approaches zero as a power 
law function of temperature in both high temperature 
and low temperature limits. 

To calculate the charge transport in this limit we in- 
troduce Keldysh technique generalized to Nambu space 

if,—G + i— Gf, + [H, G] + hoi - (2) 

where H is the Hamiltonian of electrons, including ran- 
dom impurity potentials. Icoi is the electron-phonon col- 
lision integral, (^(r, r', t, i') is Keldysh Green function 
matrix defined in Nambu space 
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P ^ R,A,K and 



0,T, 



Tz is z-component of 



Pauli matrix in Nambu space. Let Gp {e, r, t) be Keldysh 
component of the semiclassical Green function defined as 

Gp(e, r, t) = — / d^pdr'dt' exp(iet' — zp • r') 



G(r + r72, r - r'/2, t + t' /2, t - t72) (4) 

where = p'^/2m — e^^, ep is the fermi energy. In diffu- 
sion limit, Gp — G(e, r, i) + n • l\7G{e, r,t), n is the unit 
vector along p. In the adiabatic approximation, the so- 
lution of Eq.2 for Keldysh component is written in term 
of two noncquilibrium functions /i, /2 and G^'^, 

G^(e, , r) = (G^(6, r) - G%, r))(l - 2^^ (e)) 
+G^(e,r)(/2 -f r,/i) - (/s + r,/i)G^(e, r), (5) 

with np as Fermi distribution function, /i(2)(e) is an 
even(odd) function of e. 

Following Eqs.2,5, in the first order adiabatic approx- 
imation, the equations of /i^2 and G^ take forms of 

DV.i?2(.,r)V/2(.,r) = ^ Pde'^"^'^'^ 



de Jq dt 
V-i^i(e,r)V/i(e,r) = 

ze[r„ G^(e, r)] + ■ (G«(e, r)VG^(e, r)) = 0, (6) 

where D is the diffusion constant in bulk metals; v{t, r) = 
Tr(G«T, - TzG% Z?i(e,r) = Tr(l - G^r.G^r,), 
D2(e, r) = Tr{l — G^G"^), Tr represents trace in Nambu 
space. Here we neglect /I's contribution in first equation 
in the leading order oi T/ep- To derive the right hand 
side of the first equation in Eq.6, we use an identity 
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which is valid when the external parameters vary little 
over the characteristic length scale of G^'^. The second 
term in the right hand side of Eq.7 vanishes following 
the sum rule. Inside the leads, quasiparticles remain in 
equilibrium and /i(e,r) = /2(e,r) = at the lead-metal 
boundaries. At SN boundaries, G^ = 0, = exp(ix) 
at e ^ A. Electron-phonon collisions are neglected in 
Eq. 6 in the low temperature limit when the diffusion 
relaxation dominates. 

Finally, in the leading order of 1/kpl, the current can 
be expressed as 

J^Diyo f de[-^7^2(e,r)V/2(e,r)-KZ?i(e,r)V/i(e,r)], 

(8) 

with fo the density of states in bulk metals. 




Fig.l. Geometry of the sample, W d. Shaded bars 
stand for superconductors; the patterned region repre- 
sents a disordered normal metal. 

For the quasi 1 — D geometry shown in Fig.l where 
ly, L > d > ^JlTfK, G^'^(x,?/) of junction 1,2 can be 
approximated with G^-^ijj) calculated in the leading or- 
der of {d/W)'^{L is the length of the sample along x- 
direction) at given xi,2- This allows us to integrate Eq. 
6 along y-direction independently. The resultant equa- 
tions /i,/2 and G^'"^ have only x-dependences. From 
Eq.6 and its boundary conditions, we find fi — and 
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dt 



7t)-^]- (9) 



Here we introduce X{£,x) — dx' de'i'{e' ,x'), and 
p-i(e, x) = D2{e, x) dx' D^\e, x'). 

We want to emphasis that G^'^ depends on time only 
through boundary conditions Xi,2, Xi,2 G [— oo,+cx)]. 
Wave functions are periodical in such a 2D space, and 
G'^'^(Xi,X2) = G^'^(xi + 27r,X2 + 27r). Furthermore, 
time reversal symmetry of the ground state at x = 
requires that physical quantities like D12 or v are even 
functions of x- Evidently, it is more convenient to in- 
troduce a compact 2D space M2 of ga, ga — cosXq, 
a = 1,2. 5„ e [-1,1]. 

Following Eq. 6, in the adiabatic approximation, the 
noncquilibrium distribution function depends on time lo- 
cally through Xi,2(0i dtXi,2{t) or g, 5tg. Correspond- 
ingly, the time dependence of p(e,x), X{e,x) defined af- 
ter Eq. 9 can be expressed through their dependence on 
g. Thus, the time derivative of X is represented by a 
one-form defined in such a 2D parameter space. Substi- 
tuting Eq.9 back into Eq. 8 and using Stoke's theory to 
carry out the integral of the one-form along a close tra- 
jectory C in M.2 [Q, we obtain the charge transport along 
x direction per period in term of a two-form defined in 
M2 



TTapdga A dgfj 
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(10) 
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where S is the region enclosed by trajectory C in Ai2- 
In the leading order of (d/W)^, tTq/? only depends on the 
sensitivities of D2, v with respect to g, which reflects the 
degree of deformation of coherent wave packets at point 
g in M2- 



7ra/3(e,g) = 



dgn dga 

S(e,5a) = / de'i^{e',ga), 
Jo 

2/3 (e,g) = 



-02(6,57) 



7]o + (1 - 2ryo)]-' x 



^ ^ ^2[e,g'r) 

with 02,1^ calculated at given ga in the limit W ^ d. 
Vo = W/L, rjp = (a;/3 + W/2)/L, xp is the x-coordinate 
of the center of (3 region. 0{x) is a step function. In 
the absence of phase rigidities, 'E,'E. do not depend on g 
and TTai3 = 0. Q naturally represents a connection over 
a tangent fiber bundle TM , where the fiber is one- form 
Yi{d'E./ dg)dg embedded in the tangent manifold of M.2 
and the base manifold is A^2- Such a geometric point of 
view was first emphasized in |12 13] for isolated quantum 



mechanical systems and generalized to statistical systems 
in Q, where kinetic processes are important. Generally 
speaking, Q is not proportional to the area enclosed by 
trajectory C in the base manifold A^2 because tTq^ de- 
pends on g. Absence of the area law, which usually holds 
in a classical system, signifies strong correlation effects. 
Following the normalization condition we rewrite 
— cos9{e),F^ — exp(ix(e)) sin6'(e), where 9 — 
9 1 -\- i92, which is obtained as a solution to the last equa- 
tion in Eq. 6 ||. Di[e) = cosh^6'i(e), D2{e) = cos2 6'2(e) 
and v{e) — cos 6*1(6) cosh 6*2 (e). It is important to notice 
that 6*1 (e) 7^ 7r/2 or 1)2, v ^ Q only when e > Eg, where 
Eg is a function of x H] 



Eg{x<n)=E, 



^2(1 



C'iX^),X<7r; 



(12) 



Here Ec = D/d"^, Ci = 0.91, C2 = 3.122, C3 = 2.43. Fol- 
lowing Eq. 12, the density of states has a gap, which 
closes only when x = tt. At energies much larger than 
Ec, the phase dependence of ^1^2 decays exponentially as 
a function of energy. 



(13) 



9i,2 =exp(-./—)(l -1-51^2). 

V 

Taking into account these results, we find 



7ra/3(e,g) = 



'[9{ria-r/f))+r/f)-l],e:$> Ec 
F{ga, + l,gf3 + l,e/Ec),e<$:Ec. 




Fig. 2. a). Distribution of TTap in A^2 at kT ^ Ec- tTq^ 
is nonzero only in the shaded region of size kT/Ec- Tra- 
jectory a represents the case S = t:/2 and trajectory /3 
for S ~ kT/Ec. b). Q(T) as a function of (5 at fcT > Ec 
(dashed fine) and kT <C (solid line), c), d), e), f) are 
plots for {vi{T, t), V2{T, t)) at different 5 and at kT Ec- 
c). 5 > kT/Ec, d). 5 approaches kT/Ec, e). 5 ~ kT/Ec, 
f). 5^0. 

F is zero when \go,.f} + 1| S> [e/EcY and approaches 
unity when \ga,i3 -I- 1| ~ {e/EcY- Eqs. 11, 14 indicate 
that 7ra^(e,g) is exponentially small and is a simple si- 
nusoidal function of x at e 3> i5c- In the opposite limit, 
the phase dependence of tToli3 depends strongly upon the 
high harmonics. 

At high temperatures, the main contribution to tTq^ is 
from electrons of energies of order Ec ^ kT since con- 
tributions of electrons with energies of order kT are ex- 
ponentially small. Substituting Eq.l4 into Eq. 10, we 
obtain the charge transport at high temperature due to 
deformation of quasiparticles, 



Q{T) = N 



Ec 



Ec_ 
kT 



Fi{S). 



(15) 



N is the number of electrons inside the sample. Fi{6) is 
a smooth function of 6, as shown in Fig.2b(the dashed 
line); Fi{S = 0,7r) = 0. The small factor Ec/ep is from 
the asymmetry of the electron-hole spcctrums at Fermi 
energy; Ec/kT originates from the fact that only wave 
packets of this small fraction of quasi-particles can be 
deformed. When kT is larger than the Josephson cou- 
pling energy, the thermal phase slippage takes place and 

g(r) = 0. 

At low temperature, tTq,^ is determined by e ^ kT <C 
Ec- Following Eq. 14, it is nonzero only in a small region 
where \ga,i3 + 1| is smaller than kT/Ec{-€i 1), as shown 
in Fig2.a. At zero temperature, tTq,^ becomes completely 
localized at ga ~ gp = —1 and equals zero in the rest 
fl4'^^ plane. The localization of tTq^ in M2 indicates 
very different phase dependences of charge transport in 



3 



the low temperature and high temperature Umits. Con- 
sider the case 5 = it/2. The corresponding trajectory a 
in Fig. 2a does not contain the small shaded region where 
TTa/s 7^ 0. Following Eq.lO, Q equals zero in this case. To 
have nonzero charge transport, trajectory C is required 
to enclose a small portion in M.2 of radius kT/Ec where 
TTa/S 7^ 0, and 5 needs to be close to kT/Ec- 

Physically, charge can be transferred only when the 
density of states at energy kT in regions 1,2 (see Fig. 
1) is changing simultaneously. It can be demonstrated 
more explicitly by directly looking at deformation of wave 
packets, which can be characterized by the trajectory of 
the vector {vi{T,t),V2{T,t)) in a 2-D space, as shown in 
Fig. 2. c). d). e). f). Here vo^i,2{T,t) are the density of 
states of energy kT in region 1,2 at time t. According to 
Eq. 12, the energy gap Eg in the quasiparticle spectrum 
closes only at x = tt and the density of states at low en- 
ergy remains zero for most of the period unless x is in 
the vicinity of tt. When 5 ^ kT/Ec, vi changes while 
1/2 remains zero and vice versa. The trajectory there- 
fore repeats itself when the time evolves and the area 
enclosed is zero. When 5 approaches kT/Ec, V2 becomes 
nonzero while v\ is changing. The trajectory starts to 
enclose a finite area. The area reaches a maximum when 
5 ^ kT/Ec- When 5 decreases further towards zero, vi 
and V2 start to change in phase; the trajectory becomes a 
narrow strip with its width proportional to 5. As J = 0, 
the trajectory again repeats itself in the 2D space and 
the area enclosed is zero. The trajectory of (5 < and 
that of (5 > obeys mirror symmetry, i.e. the trajectory 
of 5 < can be obtained by reversing the direction of the 
trajectory oi —6. As a consequency, 

«r) = A'(^^f)V.(..f). (16, 

where 

r |(5|5o, |5| « Jo; 
F2{5,5o) = sign5 { 5l,\5\^5o; (17) 

{ exp(-|J|/5o),|<5| »<5o. 

The contribution to the charge transport in the limit 
I (5 1 ^ 5q is from the quasiparticles of energy of order 
5 X Ec, population of which is much smaller than that 
of typical excitations of energy kT. Though there is an 
energy gap in the spectrum, charge transport is not quan- 
tized in this case because kT can be higher than Eg{x) 
when X ^ 0. 

From Eqs. 15, 16, we find the charge transport van- 
ishes as the temperature goes to zero because of the 
electron-hole symmetry and decreases as 1/T at high 
temperature limit due to the suppression of electron-hole 
coherence. Q as a function of temperature develops a 
maximum at kT ^ Ec. At asymptotically low temper- 
ature, one can no longer neglect the finite dwell time of 
the electrons in region 1,2. Q{T) saturates at the value 



N{d/W)^{Ec/eF)'^ when kT - Ec{d/W)'^. Another type 
of quantum interference effect, i.e. mesoscopic fluctua- 
tions is negligible in the present case because the dimen- 
sion along z direction in Fig.l. is assumed to be infinite. 

In conclusion, we find that electrons can be transferred 
out of dots with the help of the phase rigidity of coherent 
wave packets. In the presence of small amount magnetic 
impurities, the gap of quasi particle spectrums closes as 
far as phase differences are in the vicinity of tt, the size 
of which is determined by the pair breaking rate. In this 
case, TTa/s would be delocalized in Ai2 space and the pro- 
nounced peak of charge transport as a function of 6 as 
shown in Fig. 2b will be smeared out. When the depair- 
ing rate is further increased, coherence effects become 
exponentially small and Q{T) w 0. 

Finally, in the metallic limit studied in this letter, 
quantum phase slippage effects are negligible. 
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